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SYMBOLS 


,a^ numerical constants for parabolic velocity profile 

,bj numerical constants for parabolic velocity profile 

phase velocity of vortex 
amplitude of vortex 
circulation constant 
order of magnitude 
inner region independent variable 
Reynolds number of undisturbed flow 
time 

velocity in streamwise direction 
free-stream velocity 
velocity in transverse direction 
streamwise coordinate direction 
streamwise coordinate in convective frame 
streamwise coordinate in inner region 
transverse coordinate 

transverse coordinate in convective frame 
transverse coordinate in inner region 
circulation of vortex 
boundary layer thickness 
small matching parameter 

small parameter governing size of inner inviscid region 
kinematic viscosity 
arbitrary phase angle 
stream function 


vorticity 



Subscripts ; 


c 

i 

np 

o 

P 

s 

0,1 ,2 


center of vortex 
inner region 
nonparallel flow 
outer region 
parallel flow 

+ when analysis is for y > y ; - when analysis is for y < y„ 

c c 

order of approximation 


iv 



SUMMARY 


Tile flow field of a vortex in a viscous shear flow is found by constructing a 
uniformly valid asymptotic expansion consisting of an inner solution field repre- 
sented, to lowest order, by a two-dimensional, nonlinear, inviscid Stuart vortex and 
an outer solution field represented, to lowest order, by either a two-dimensional 
parallel or self -similar viscous flow. The technique involves scaling both the 
transverse and streamwise coordinates in the vicinity of the vortex as well as allow- 
ing for a "slow" variation of the outer viscous flow. Criteria are established for 
both the size of the vortical structure and proximity to the boundary surfaces. Tiie 
composite solution is a consistent mathematical picture of the flow field at a fixed 
streamwise location as the vortical structure evolves past this point. Such a formu- 
lation is also useful in the specification of boundary or initial conditions in 
numerical fluid dynamic calculations, where an inconsistent setting of these condi- 
tions leads to spurious results for rather long computation times. 
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1 . INTRODUCTION 


Of general interest in a broad class of wall bounded flows is the dynamic evolu- 
tion of vortical structures through the flow. In general, these structures are 
three-dimensional and their downstream evolution is exceedingly complex. In fact, as 
yet no overall mathematical description of such entities has been formulated. It is 
one of the intentions here, however, to establish a framework, based on first princi- 
ples, which may form a basis for more detailed analytical studies. Another motiva- 
tion concerns the establishment of boundary and initial conditions in numerical 
experiments. In such problems it is of great importance that any distribution of 
values, for the dependent variables, be consistent so that no extraneous or nonphysi- 
cal results contaminate the solution field. Ihis is especially true when such a 
variable as pressure is required from the numerical solution. 

The mathematical framework that is used here is the method of matched asymptotic 
expansions. An inner solution field is constructed which consists of a two- 
dimensional vortical structure. The outer solution field is taken to be an otherwise 
undisturbed laminar two-dimensional parallel or self -similar viscous flow field. In 
order to allow for a finite-amplitude disturbance field, both the transverse and 
streamwise coordinate directions in the inner region are scaled by a suitably defined 
small parameter. Since such a short wavelength solution field must, in general, be 
coupled to the more slowly varying outer flow, a two-variable ejqiansion procedure in 
the streamwise direction is necessitated. Ihe coordinate basis allows for the com- 
plete description of the inner flow field. In the outer field, the scaled transverse 
coordinate is replaced by a transverse coordinate characterizing the overall size of 
the boundary layer flow, and the streamwise coordinates are left unaltered. 

It has been argued that the dominant dynamics in the vicinity of the vortex are 
inviscid (e.g., ref. 1) and that viscous effects, which are omnipresent, play a sec- 
ondary role in the overall features of the vortex. With this as a basis, an exact 
solution to the appropriately scaled nonlinear, inviscid vorticity transport equa- 
tion, the Stuart vortex solution, is taken as the inner solution field to lowest 
order. It is found that, to lowest order, this solution is easily matched to a class 
of parallel and nonparallel viscous dominated flows, Ihe composite solution field is 
constructed from the lowest order inner and outer solutions, and higher order equa- 
tions are used to identify the inner variable scaling parameter. Finally, restric- 
tions on the proximity of the vortex to boundary surfaces are derived from the 
requirement of zero component velocities on the boundaries. 


2. INNER SOLUTION FIELD FOR EMBEDDED VORTEX 

Consider a two-dimensional vortical structure evolving in an otherwise undis- 
turbed two-dimensional viscous flow. Ihe underlying dynamics of the vortex is invis- 
cid and its evolution is nonlinear. With these basic requirements, one can deduce 
from first principles a vortical structure which is embedded in an otherwise viscous 
dominated flow. As indicated in section 1, the method of matched asymptotic expan- 
sions is used to construct a uniformly valid distribution of flow variables in a 
class of parallel and nonparallel wall bounded laminar flows. 

An appropriate frame of reference must first be established in the flow. At 
some time t^, the vortex with center at (x^,,y^) is assumed to be moving uniformly 



with constant phase velocity c. The value c used is taken as the Eulerian veloc- 
ity at the point (x^,y^) in the undisturbed flow field. The vortex flow is then 
stationary in the coordinate system 


x' = X - ct (1 ) 

y' = y - Yc (2) 


where x' and y' are the streamwise and normal coordinate directions, respec- 
tively, and x' = 0 • when x = x^ and t = t^. Such a description allows for the 
constructed vortex to be placed in a steady viscous outer flow field which is either 
independent of x itself or where the flow has reached some self-similar state. In 
addition, the coordinate lengths x' and y' , appearing in equations (1) and (2), 
have been scaled by some characteristic length of the flow such as channel height or 
local boundary layer thickness. This, of course, is necessitated by the need for 
0(1) scaling in both the inner and outer solution fields. The physical problem of 
interest requires vortical structures of small spatial extent relative to the charac- 
teristic scales of the overall flow. Such a requirement dictates a rescaling of the 
independent variables in the region dominated by the vortical structure. This yields 
a stretching of the y' coordinate and the decoupling of the streamwise dependence 
into a slow variable, representing the slow divergence of the overall flow, and a 
fast variable, representing the locally rapid streamwise variation of the disturbed 
vortical flow. The governing differential equations for the embedded vortex motion 
are then written in terms of the inner variables: 


X 


$ 




(3) 


where y is a small scaling parameter which is a function of the inverse Reynolds 
number and is chosen such that in the vicinity of the vortex, x and y are 0(1). 
In general, x can be expressed as a series expansion of x' , which would account 
for any wave-number variations at higher order; however, to the order of the analy- 
sis, these higher order terms can be neglected. In order to apply equation (3) to 
slowly diverging flows and maintain a uniform phase velocity along the symmetry line 
y = 0, it is necessary to assume a dependence of y^ on the slow variable x' , 


y = ^[y - 


(4) 


This two-dimensional problem is formulated in terms of the stream function and 
vorticity, and in terms of the previously defined inner variables, the governing 
differential equations can be written as 


u 


1 ii 
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+ -=~ 

3x2 gp 

and R is a Reynolds number formed from a characteristic length and velocity scale 
of the flow, \|; = i|;(x',x,y) is the stream function, u = u(x',x,y) and 
V = v(x',x,y) are the streamwise and normal velocity components, respectively, and 
Q = £2(x',x,y) is the vorticity. Coupled with the governing differential equations, 
an appropriate asymptotic expansion of in the inner vortex region is required. A 

suitably general expansion for i|) is given by 


\|» = \pi(x',x,y) = uR“^ ,x,y) + vi2j^-1 (x' ,x,y) 


+ vj2r-1 \(;^ 2 (x' , x,y) + 0(y'*) 


(9) 


where the usual requirements for a valid asymptotic sequence are assumed to hold 
Before obtaining the governing differential equations at the various orders of 
approximation, it is necessary to quantify the factor appearing in the equations 
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In anticipation of the Blasius boundary layer example to be shown in section 5, 
one can estimate the size of dy^/dx' for this boundary layer flow. If the trans- 
verse coordinate is scaled by the boundary layer thickness at some downstream loca- 
tion Xq, that is. 


6 


0 



( 10 ) 


where U is the free-stream velocity outside the boundary layer, then y^ is anal- 
ogous to the usual Blasius similarity variable with 0 < y^ < 1 . Thus the factor 
dy^/dx' can be written, to lowest order, as 


dy 

dx' 


5 

2 


y R 
^ c X 



(11 ) 


where = — — in this boundary layer case. With this specification completed, 
the governing differential equations at the various orders can be derived. 

The 0(u) equation is obtained by substituting the expansion equation (9) into 
equations (5) through (8): 


■f^('l'io - cy) 

3y 


9^i0 

9x 


9^i0 
9x 9y 




- 1/2 


9y 


— 


yR 


iO 


( 12 ) 


where 


= -V il;.„ 
xO ^xO 


(13) 


Note that in the present formulation, the lowest order vorticity is not simply deter- 
mined by gradients of the streamwise velocity but is determined by gradients of both 
the streamwise and transverse velocities. In order to be consistent with the initial 
physical assumptions concerning the dominant inviscid behavior of the vortex motion, 
it is appropriate that the limit (pR)“^ -»■ 0 be chosen in equation (12). The 
0(y) equation reduces to 


9^i0 ^^iO 9^i0 ^^iO _ 

9y 9x 9x 9y 
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or 


3 (4^io ) 

3 (y,x) 


( 15 ) 


where 


^iO = 'f'iO - cy, JJiO = 

The set of equations, equations (13), (15), and (16), were originally studied by 

Stuart (ref. 2) in the context of finite-amplitude oscillations in mixing layers. 

The solution to the set, appropriate to the physical problem at hand, is 

’l'io(^»y) = cy + K ln[cosh y + Cq cos(x + (Ji)] (17) 


S)i0(x,y) = J^io(^io) = -k[i - '^o^] (18) 

where c^ is a constant constrained by (ref. 1) 

o < c <1 (19) 

o 

and K is a constant associated with the circulation of the vortex having bounding 
contour, -<j) < x < 2ir - <j) , -«> < y < <», and (f) is an arbitrary phase angle. The 
circulation of the vortex, as given by equation (18), is easily determined (ref. 2) 
as 


r = -4 ttk 


(20) 


The distribution given in equations (17) and (18) is rather general, in the 
sense that the constant c^ governs the gross features of the vortex in a rang;e 
varying between a point vortex distribution (c^ 1 ) and an unperturbed hyperbolic 

tangent velocity profile (c^ = 0) . The stability characteristics of this vortex 

distribution have been studied in reference 3; however, the intent here is not to 
conduct such a stability study but to examine the composite flow structure at some 
point in space and time when such a vortex is evolving through the flow. This is 
also distinct from recent work by Pierrehumbert and Widnall (ref. 4) and the earlier 
work of Browand and Weidman (ref. 1), where such vortical distributions were studied 
in connection with free shear layers. There the objective was to examine the inter- 
action dynamics and critical parameters of interacting vortices which spanned the 
width of the shear layer . 
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This lowest order solution leaves the two parameters y and K unspecified. 
Higher order inner solutions are needed to fix y, and the matching process will help 
structure the outer solution field. The O(y^) differential equation is 


3(iJ;i1,nio) 9(^i0»^il) 

:: — ~ + :: — z 

3(y,x) 3(y,x) 


1 


/_ 9(?i0'^i0»xy) 3J2iol 

yRl /2 
X 

-■ 2j 

S2 » _ + c » / 

3(y,x) 3y J 


+ 


1 

y^R 


v2nio 


(21 ) 


where 


J2i1 = , 


^iO,xy 


3x 3y 


( 22 ) 


The question that arises at this point is the choice of the appropriate limit for 
(yR^'^)~^ . Since the differential operator on the left-hand' side of equation (21) 
is different than the lowest order differential operator, it is consistent to solve 
the homogeneous problem at this order; therefore. 




0 , 


y » R 


- 1/2 


(23) 


is the chosen limit. The homogeneous form of equation (21 ) is readily solved by 
introducing the new independent variable 


Q = 



/K 


(24) 


and transforming to the new coordinate system (Q,x). 


The resulting equation is 


■^[F] = 0 (25) 

3x 


where 


F(Q) = n., - 2(1 - c ^)Q^t-. 

i1 o i1 


( 26 ) 
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and, in general, 0^^ and can both retain some periodic dependence on x. It 

is necessary to continue to the next order to further identify £2. and and 

uniquely determine y . 

3 

The o(y ) equation can be written as 


3 (ipiO f^^i2 ) 9('l'i2'^i0^ 

3(y,x) 3(y,x) 


3 ,£2^1 ) 1 r ^ ~\f 3(tio>?iO,xy) 

3(y,x) y2j^1/2|J 2 3(y,x) 


+ 



+ 


1 

y^R 


V2£2. 


iO 


(27) 


where 


i2 



(28) 


Once again transforming to (Q,x) coordinates would place equation (27) in a more 
tractable form. In light of the O(y^) equation (eq. (25)), the right-hand side of 
equation (27) cannot retain any secular terms, that is, any periodic functions of 
X. In the absence of any decision on the limiting values of the terms continuing y 
and the Reynolds number, it is necessary to examine the functional behavior of 


3 ('I'ii f^il ) 
3 (y,x) 


or 


3 (’I'ii ) 
3 (Q,x) 


G(Q,x) 


(29) 


in (Q,x) coordinates. If = ’I'il (Q/X) and £2^1 = £2^1 (Q,x), then equation (29) 

becomes 


liilf 

3x L 


+ 4Q(1 



(30) 
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where ' denotes differentiation with respect to Q. In order for this term to be 
independent of x it is necessary for , and therefore , to be functions only 

of Q. The structure of these two functions can be further delineated by referring 
back to equation (26). Equation (26) can be expanded to give 


F(Q) = q(i - ) + Q'l'il - - f(Q.’l'iT) 

+ 2 Cg cos(x + (})) f(Q^il;^^) 


(31 ) 


where 


f(Q,i^.^) = Q^r[^ + (32) 

Clearly f(Q,\j;. ) must equal zero and the resulting solution is = Constant. 

Without any loss of generality in the inner solution field, it is sufficient to take 




i1 


= n., = 0 


(33) 


Now the limit processes in equation (27) can be performed. Since inhomogeneous terms 
should appear at this order, it is only necessary to determine whether nonparallel or 
purely viscous effects predominate. In the absence of any nonparallel effects the 
viscous limit dominates, that is. 


y 


p 


R 


-1/3 


If the outer flow is slowly diverging, the appropriate limit is 


(34) 


np 


= y_ = R 


-1/3 


(35) 


since 


1/2 

R > MR 

X 

Even though this analysis was not motivated by questions of stability, it is worth 
noting the viscous layer size obtained here (eq. (34)) and to compare these neces- 
sarily nonlinear results with previous nonlinear stability studies, such as those 
motivated by the ideas of Benney and Bergeron (ref. 5). In those studies nonlinear 
effects dominated in the critical layer and viscous effects were secondary. This, of 
course, is consistent with the structural ideas presented in this work. Finally, it 
should be pointed out that the nonparallel effects are sustained by terms originating 
in the nonlinear advection terms. In the absence of such stress gradient effects. 
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the slow variation of the outer flow would not be felt by the vortex to this order of 
approximation . 


3. INTERMEDIATE MATCHING REGION 

In section 2, the inner solution field was established. This solution must be 
formally matched to an appropriate outer solution in order to form a uniformly valid 
composite expansion. The process involves the introduction of suitable intermediate 
variables (ref. 6) given by 


y 


e 


y - y^(x’) 
e 


(36) 


where 


e 

— ^ oo as e ^ 0 

y 


and 


y = 


eye 


as 


y - yc| ^ 0 


since the vortex is embedded in the flow there are two distinct matching regions: 
the region above the vortex y ->■ » and the region below the vortex y ->■ -<». If the 
lowest order inner solution is written in terms of the intermediate variables, the 
inner solution (eq. (17)) becomes 


If 




'I'i = cygE + pK In l^cosh — — + Cq cos(x + (j))j + O(p^) 


<*>, the two matching limits of equation (37) can be written as 


(37) 


1 );^^ -»• (c + K)Eyg - pK In 2 + pK In j^1 + c^e ^^^^^^/cos(x + 4>)j 
+ O(p^) -»■ 


(38) 
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( 39 ) 


-»• (c - K)eyg - yK In 2 + yK ln|j + CqS + <|)) 

3 \ 

+ o(y ) 1-^ -»j 

These stream function values will indeed match to the o( 1 ) outer stream function 
distribution, and the contributions to the respective outer flow regions are 


1 

ll+“ 

(c + K)ey^ 

(40) 

1 -»■ 
i| -« 

(c - K)ey^ 

(41 ) 


It is important to note that in the present formulation, the periodic streamwise 
dependence is lost in the matching process because of the transcendentally small 
effect. This eases the restrictions on the interpretation of the present formulation, 
since the streamwise range of the vortex disturbance is limited here to a single 
wavelength. In addition, this allows the inner solution to be matched locally, point 
by point, to an outer flow which is either parallel or self-similar, but nevertheless 
slowly diverging. In the next two sections, an example of each type of flow will be 
given. 


4. OUTER VISCOUS FLOW FOR PARALLEL FLOW CASE 

The independent variables appropriate for the general outer problem are x' , x, 
and y and are assumed to be scaled to 0(1). For the moment, the derivations will 
be rather general in order to lay a proper framework for both the parallel flow case 
considered in this section and the nonparallel case to be considered in the next 
section. The general form of the equations for the two-dimensional outer flow prob- 
lem can be written as 


u 


3y 


(42) 


1 

3x' y a“ 


(43) 




iil _ ill 

9y2 9 x' 2 p2 9x2 


2 92i|; 

^ 9x 9x' 


(44) 


9(2 

+ 1 M 

. 11- + 

1 11 

lil . 

- 1 llR + 

92(2 

^1 92(2 ^ 

2 92(2 

9x' 

M 9x_ ■ 

9x' ^ 

V 9x 

9y ■ 

R 9y2 

9x' 2 

y2 3x2 

^ 9x 9x'_ 


( 45 ) 
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Recalling the form of the outer limit of the inner solution, it is seen ;^at the type 
of outer flow that is considered can be explicitly independent of both x and x' . 
Then, if the outer flow is parallel or even dependent on some similarity variable, 
the lowest order outer solution and inner solution can be matched. The example prob- 
lem to be considered in this section is plane Poiseuille flow. The outer dependent 

variables are independent of x and x' to lowest order, that is 

= Tp (y) (46) 

o os 

where s is + for y > y^, and s is - for y < y^,. Equation (45) reduces to 

the trivial form 




os 


= 0 


3y 


(47) 


with solutions 




o+ 




1 ) ^ + 


1) + a3 


(48) 







(49) 


where the variables have been scaled by the channel height h and the maximum undis- 
turbed inflow velocity U . The boundary conditions imposed at lowest order are 

00 


(y = y„) = (y = y) = o 

0 + c o— c 


(50) 


os 


2(1 - y ^ ) - |-( 1 
c 3 


. o 2 4 3 

i-2yc -^3^0 


^c ’ 


( s + at y = 1 ) 


(s - at y = 0) 


(Undisturbed values) (51) 


3ip 


o+ 


dtp 


o- 


3y 


y=1 


3y 


= 0 


(52) 


y=0 
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In anticipation of the matching, the value for the vortex phase velocity needs to be 
determined. This velocity is the velocity at y = y^, if the flow did not contain 
the vortex. In the channel flow case 


c = 4y^(1 - y^) 


(53) 


If the boundary conditions are applied to equations (48) and (49), the constant coef- 
ficients become 


6a, 


3a, 


^ = 


(1 - y ) (1 - y ) 

c c 



(54) 


a.^ = 0 . 0 , 


= 0.0 


(55) 


a 


3 





(56) 


where the remaining coefficients, and bg , are determined from the matching. 

The outer variables (1 - y) and y can be written in terms of the matching variable 
y^ and are given by 


1 - y = (1 _ y^) _ Sy^ 


(57) 


y = 


^c 




(58) 


Substituting these forms into equations (48) and (49), one can obtain and 

in terms of powers of ^y^* From section 3, it is seen that the matching occurs with 
coefficients of ^y^* This process yields the coefficient values 


6a 


^0 


3 ^ 2(c + K) 


(1 - y ) 

(1 - y ) -^c 

c 


^^3 2(c - K) 

^0 = - y 


(59) 


(60) 


The remaining coefficients a^ and b^ are then determined from these values of 
aQ and b^. Note that the coefficients are all bounded in the present formulation. 
This results from the fact that the inner region is assumed to be of 0(y) in size; 
therefore, the vortex should be at least a distance of 0(y) away from the bound- 
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aries to be consistent with the basic assumption about the inner region. A uniformly 
valid composite solution, to lowest order in the dependent variables, can be formed 
by summing the lowest order inner and outer solutions and subtracting out the common 
part (ref. 6). The resulting composite solutions are 


■r^d-y)^ - (1-y )^] - -^[(l-y - (1-y )^] - K(y-y ) + KR ln[cosh y + c cos(x+4i)] (y>y ) 

2 OoOOO o c 


'I' “ < 


> (61 ) 


•^^y^-y^^] + ^ K(y-y^) + KR ^ ln[cosh y + cos(x+it') 


(y<y^) 


r 


o ~ J '-O 

-agd - y) + —(1 - y)^ - K + K tanh ySI + 


{ 


c„ cos(x + 


cosh y 


(y > yc) 


u “ N 


> (62) 




b -1 ~ J '^o cos(x + <{)) 

boy + -7^ y2 + K + K tanh y 1 1 + = 

2 cosh y 


<y < Yc^ 


V « K 


Cn sin(x + (()) 
cosh y 


{ Cq cos(x + <j))| 1 
+ -= / 
cosh y J 


(0 < y < 1 ) (63) 


n « -KR^^^[l - Cp^] sech^y^l + 


^ T fa,- - a (1 - y) 
Cq COs(x + (|))|''^ ' ® ' 


cosh y 


bo + b^y 


(y > Yq) 

(y < y^) 


(64) 


where x, y, and (f) were defined in section 2, u = R ^ l^as been used, and the 
constant coefficients were defined earlier in this section. It should be pointed out 
that the vorticity is discontinuous at y = at first order, 0(1); however, to 

zeroth order, 0(R^^^), the vorticity is continuous across the channel. If higher 
order terms are included in the expansion procedure, this singularity would be 
removed. As can be seen from equations (62) and (63), the velocity boundary condi- 
tions at the wall will, in general, not be satisfied to the order of the approxima- 
tion for all values of y^ and K. The task then is to delimit the range over which 
y^ and K may vary. At y = 0 and y = 1 the velocity fields reduce to 


u(y=0) = K<1 - 


v t, 1/3 
tanh R yQ 


1 + 


Cq cos(x + (|)) 
cosh R^/^yc 


-1 


} + 0(R~^^^) 


(65) 
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v(y=0) 


Cq sin(x + (j)) 
cosh 


+ 


Cq cos(x + <{>) 
cosh Yq 


-1 


+ o(r"^/^) 


( 66 ) 


u(y=1) = -K<1 - tanh R^'^^d - y^) 


1 + 


Cq cos(x + (|>) 


cosh R^/^(1 - Yq) 


-1 


+ o(r"^/^) 


(67) 


v(y=1 ) = K 


Cq sin(x + (|)) 


cosh R^ /3 ( 1 - Yq ) 


1 + 


Cq cos ( X + ({) ) 


cosh R^ /3 ( 1 - Yq ) 


-1 


+ 0(R ^/^) 


( 68 ) 


where the accuracy of the approximation 0(R~^/^) in the velocity field is based on 
the order of approximation obtained in the inner solution field. Since the periodic 
functions are bounded, the restrictions on y^ must come from the hyperbolic func- 
tions. Note that the above does not put any restrictions on the value of K; how- 
ever, restricting K to values less than 1 would improve the accuracy with which the 
boundary conditions are satisfied. A cursory examination of equations (65) 
through (68) reveals that the hyperbolic functions are bounded as 


tanh = 1.0 - 0(R ^^^ ) 

c 


(69) 


cosh R^'^^y = 0(c R^^^ ) 

c o 


(70) 


tanh R^/^d - y ) = 1 .0 - o(R ^^^) 
c 


(71 ) 


cosh R^'^^d - y ) = 0(c R^^^) 
c o 


(72) 


where the constant c^ is included in the restrictions on the hyperbolic cosine 
factor for generality. If an 0(1) constant of proportionality is assumed in these 
order estimates and finite- amplitude disturbances are considered, then the following 
limits on can be extracted from equations (69) through (72) 


ln(2c R^^^) < y < 1 - ln(2c R^'^^) 

o « •'c » o 


( 73 ) 
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where 


2c > 1 (74) 

o 

With this determination of the limits on as much information as possible has 

been extracted from the solution to this order. The only parameter in the problem 
that has not been delineated in any concise way is the circulation constant K. 
Constraints on K may arise at higher orders, either through the matching process or 
’ the application of the boundary conditions. Nevertheless, for the present study, K 

is a free parameter of 0(1), whose value will be chosen to best suit the practical 
» problem at hand . 

As a test case, computations were carried out at a Reynolds number, based on 
channel height and maximum undisturbed inflow velocity, of 1 0 . The amplitude of the 
Stuart vortex c^ was set at 0.25 (refs. 1,3, and 4) and K was chosen to be 0.20. 
The value of y corresponding to this Reynolds number is 0.046. The results of the 
streamwise velocity distribution at various values of (x + (()) are shown in figure 1 . 
The results reflect, as (x + <()) varies from 0 to 2ir, the effect of a passing vortex 
on an otherwise fully developed parallel flow at a fixed point in space. It is 
important to note that the vortex, which is of limited spatial extent in the present 
example, affects a rather large portion of the channel flow. This result indicates 
that quantities such as surface pressure and wall shear stress are modified by the 
passage of such a vortical distribution. 

in the next section, the more relevant case of nonparallel flow is considered; 
however, the guidelines established in this section can be applied directly to this 
more general case . 


5 . OUTER VISCOUS FLOW FOR NONPARALLEL FLOW CASE 

As indicated earlier, to the order of the approximation, nonparallel self- 
similar outer flows can be examined within the present context. Consider now, as the 
nonparallel flow example, a Blasius boundary layer flow. Since the inner Stuart 
vortex solution was found to be independent of x' , the procedure used to obtain a 
composite solution for the channel flow can be applied to the boundary layer flow 
over a continuous range of x' , with the only change at each x' location being the 
location of y^. It should be recalled at this point that the nonparallel effects 
which appear in the inner region equations as dy^/dx' were approximated to lowest 
order by the constant term -5y^,R~ ' '^/2; thus, any x' variation of dy^/dx' was 
removed. The lowest order equation in the outer region is obtained from equation (45) 
by using the following approximation to ij): 

(x* ,y) (75 ) 

o os 
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where, as before, no outer flow dependence on x is assumed. The lowest order vor- 
ticity transport equation is 


os 

3(y,x') 




OS 


3x 


f2 




OS 


ay" 


(76) 


Now, if the usual boundary layer approximations are made, equation (76) can be 
written in the familiar fom 


3u 3u 1 3^u 

3x' ^ ^ 3y R - 2 
3y 


(77) 


f, 


♦ 


where R is a Reynolds number formed from the free-stream velocity and the 

boundary layer thickness at x = x' + Xq, that is 


6(x) 


vx. 


(78) 


In a completely analogous manner with the usual Blasius method, the resulting stream 
function equation is 


2il;’" + 25\b" ij; = 0 

os os OS 


(79) 


where the differentiation can be interpreted to be with respect to y in light of 
the nondimensionalization by 6(x). Equation (79) can be solved given the appro- 
priate boundary conditions. As was the case in plane Poiseuille flow, the outer flow 
boundary layer is divided into two regions, y < y^, and y > y^,. First the region 
y < y^ is considered. 

Since the solution to the ordinary differential equation (eq. (79)) will be 
determined by numerical integration, it is necessary, for matching purposes, to 
expand the outer stream function solution in a Taylor series expansion about y^, 


l|) (y>y ) 

o- c 


» t); (y=y ) + (y - y ) 
o- c c 



y=Y 


c 


+ 


(80) 


Note that the question of continuity of 3\p^ /3y across y^ does not enter into the 
present context, since the formulation of a composite solution valid over the whole 
domain is a superposition of two distinct outer fields. Each outer solution field is 
viewed as a distinct problem with appropriate boundary conditions applied at y = 0 
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and y = for y < and at y = y^ 
x' the part of the outer solution which 

’I'* (y=y^) = c - K 

o— c 

Emboldened by the success in the channel 
the perturbed outer solution is followed 
tions associated with undisturbed stream 
are 


and y ->• oo for y > y^ , Thus , at a fixed 
matches to the inner solution is 

(81 ) 

flow, the same procedure used to determine 
in the present case. The boundary condi- 
function value and zero velocity at y = 0 


(y=0) = 0 (Undisturbed value) (82) 

(y=0) = 0 (83) 

o- 


Equations (81) through (83) are sufficient to determine the stream function, as well 
as its associated derivatives, from y = 0 to y = y^,. In this example, it should 
be noted that the stream function value at y = y^, will not be 0, although it will 
be continuous since ip (y = y ) will )ae equated to = y^)* Now the region 

y > y^ will be considered. 

From the outer limit of the inner solution (eq. (79)) and the series expansion 
( eq. (80)), the inner limit of the outer solution which matches with the inner solu- 
tion is 


(y=y ) = c + K 
o+ c 


(84) 


Since the stream function is assumed continuous across y = y^,. 


, (y=y ) = (y=y ) 

o+ c o- c 

the only remaining boundary condition to be applied is at the free stream. This 
condition is given by 


lb' -)- 

^o+ 


1 . 0 , 


y -»• oo 


( 86 ) 


This, of course, is the usual free-stream condition, and the exact location where 
equation (86) holds is left intentionally arbitrary to account for any boundary layer 
thickening due to the presence of the Vortex. This extra degree of freedom was not 
present in the channel flow problem due to the rigid bounding surfaces. Equa- 
tion (79), coupled with the boundary conditions appropriate to the region of interest 
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were solved with a fourth-order variable interval Runge-Kutta method. The composite 
solution for the dependent variables can be formally written as 


r 


(y) - K(y - y ) + KR ^ ln[cosh y + c cos(x + <t>)] (y > y ) ] 

o+ c o 1 


'I' “ < 


Tp (y) + K(y - y ) + ln[cosh y + c cos(x + <|))] (y < y ) 


> (87) 


c > 


^ , ~ f ( X + <|) '' 

1^0+ (y) - K + K tanh ySl + 




cosh y 


(y > yc^ 


U « S 




+ K + K tanh y<1 + 


Cq cos(x + 4>)| ^ 
cosh y J 


(y < yc> 




( 88 ) 




Cq sin(x + ( j >) f Cq cos(x + ^ 

K = — — Si + 


cosh y 


cosh y 


V » 


Cq sin(x + <l>) j Cq cos(x + (|)) 

K Si + 


1-1 


V. 


cosh y 


cosh y 




(y > yc) 


(y < yc) 


(89) 


= -R^^^k[i - Cq^] sech'^y^l 


2 - 


=Q cos(x + 4>)\-2 fr6+(y) (y > yc) 


cosh y 




(y) (y < Yc'> 


(90) 


where x, y, and 4> were defined in section 2, y = in this case, and terms 

of 0 (r”^^^) were neglected. Note that once again the vorticity is discontinuous 
at 0(1). The restrictions on the location of the vortex can once again be derived 
from the hyperbolic functions. In light of the arbitrary location of the top free- 
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stream boundary, is only limited here to its proximity to the boundary at 

y = 0. The restrictions are 


y > ln(2 c <2c R^^^ >1) (91) 

c o o 

Once again the parameter K is free to be chosen for the particular physical problem 
at hand. 

As a test case, a computation was carried out at a Reynolds number, based on 
local boundary layer thickness and undisturbed free-stream velocity, of 10^. The 
amplitude of the disturbance c^ was set at 0.25 and the constant K was set equal 
to 0.2. For the chosen Reynolds number, the small parameter y was fixed at 0.1. 
Figure 2 shows the distortion of the streamwise velocity profile as the vortex passes 
by a fixed point streamwise location x^. Once again the effect of the vortex is 
felt throughout most of the boundary layer, even though its spatial extent is bounded 
by )i . It is also worth noting that this type of disturbance can be present near the 
top of the boundary layer. The existence of such finite-amplitude disturbances near 
the top of the layer is a requisite for the receptivity of the boundary layer to any 
free-stream disturbances. 


6. DISCUSSION 

As was stated in section 1 , the intent was to construct a uniformly valid flow 
field, to some order of approximation, that would mimic the structure of a real flow 
in the presence of a finite-amplitude vortical disturbance. The basic structure for 
the vortical field was taken to be a Stuart vortex. Such a field is derivable from 
the general vorticity transport equation when both the transverse and streamwise 
coordinates are scaled by the same small parameter. As was seen in the two test 
cases, however, the influence of the vortex spreads throughout most of the flow even 
though its dynamic features dominate a small spatial region of 0(y). In addition, 
the analysis that was presented necessarily restricts the results to large wave- 
number disturbances. 

It is beneficial at this point to place the present work in context with pre- 
vious work on vortices in the presence of solid boundaries. For example. Walker 
(ref. 7) studied the boundary layer produced by the presence of a single rectilinear 
vortex filament above an infinite plane wall. It was concluded that no steady flow 
field exists and that a recirculating eddy, which is formed in the boundary layer 
flow, causes rather explosive growth of the boundary layer itself. The problem pre- 
sented in the present paper is fundamentally different. It is assumed here that the 
viscous flow already exists and that a vortex, which is necessarily of small spatial 
extent, is embedded within this flow. In the absence of the base undisturbed viscous 
flow, the present formulation is ill-posed since it would be impossible for the 
boundary condition at the boundaries to be satisfied, in general. For these reasons, 
it is felt that the present problem closely resembles a real flow situation of a 
vortical structure evolving through an otherwise undisturbed flow. It does remain an 
open question as to the behavior of the vortex as it evolves in time and space, but 
the present work shows that a stationary situation can exist and can be described in 
a consistent mathematical framework. This framework yielded some interesting fea- 
tures which are worth noting. 
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In the analysis, nonparallel effects were accounted for by requiring that the 
entire vortex move uniformly through the flow. Thus, for the laminar boundary layer 
flow, the vortex is inclined with respect to the wall, and as was mentioned in 
section 5, this region is sustained by nonlinear stress effects. The final comment 
concerns the effect of the vortex on the boundary surfaces. Recall that the wall 
shear stress, that is 3u/3y, was allowed to vary in the calculation; only the stream 
function and velocity at the walls were taken as either undisturbed flow values or 
zero. This allowed for the effect of the passing vortex to be felt in both the wall 
shear stress and pressure signatures . 

The work presented in this paper is intended as a first step in a more complete 
and accurate mathematical description of vortical structures, three-dimensional as 
well as two-dimensional, embedded in a broad class of viscous outer flows. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
October 31, 1983 
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